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ABSTRACT 

In the neutron-rich internal shocks model for Gamma-ray Burts (GRBs) , the Lorentz factors (LFs) of 
ions shells are variable, so are the LFs of accompanying neutron shells. For slow neutron shells with a 
typical LF ~ tens, the typical (3— decay radius reads Rp, s ~ several x 10 14 cm, which is much larger than 
the typical internal shocks radius ~ 10 13 cm, so their impact on the internal shocks may be unimportant. 
However, as GRBs last long enough (T 90 > 20(1 + z)s), one earlier but slower ejected neutron shell will 
be swept successively by later ejected ion shells in the range ~ 10 13 cm — 10 15 cm, where slow neutrons 
have decayed significantly. We show in this work that ion shells interacting with the (3— decay products 
of slow neutron shells can power a ultraviolet (UV) flash bright to 12th magnitude during the prompt 
7— ray emission phase or slightly delayed, which can be detected by the upcoming Satellite SWIFT in 
the near future. 

Subject headings: gamma rays:bursts — radiation mechanisms: nonthermal — shock waves 



1. INTRODUCTION 

As realized by many authors, the fireball of Gamma- 
ray bursts (GRBs) may contain significant neutron com- 
ponent in essentially all progenitor scenarios (e.g., Deri- 
shev, Kocharovsky & Kocharovsky 1999a, b; Beloborodov 
2003a). The dynamics as well as the observable signa- 
tures (the neutrino emission from the proton-neutron de- 
coupling and the changes in GRB lightcurves due to neu- 
trons preceding the ion shell or catching up it) of a rel- 
ativistic neutron-rich fireball have been investigated by 
Dcrishev, Kocharovsky & Kocharovsky (1999a, b) firstly. 
Since then, many authors have had their attention on the 
neutron- fed GRBs. If the neutron abundance is compara- 
ble to that of proton, the inelastic collision between dif- 
ferentially streaming protons and neutrons in the fireball 
will provide us observable 5-10 GeV neutrinos (Bachall & 
Meszaros 2000, hereafter BM00; Meszaros & Rees 2000). 
The baryon-loading problem in GRB models can be ame- 
liorated if significant fraction of baryons confined in the 
fireball are converted to neutrons (Fuller, Pruet & Abaza- 
jian 2000). Further investigations on the observation im- 
plications have been presented in Pruet & Dalai (2002, 
hereafter PD02) and Beloborodov (2003b, hereafter B03). 
In PD02, the short GRBs are assumed to be powered by 
external shocks and the neutron component lags behind 
the ion one. As the ion ejecta has been decelerated sig- 
nificantly by the external medium, the decayed neutron 
ejecta catches up the decelerated ion ejecta and powers 
keV emission. In B03, the neutron-ejecta keeps ahead the 
ion-ejecta. The f3— decay (n — > p + e~ + f c ) products share 
their momentum with the external medium immediately. 
As a result, the external medium have been accelerated to 
a ultra-relativistic velocity. The interaction between the 
ion ejecta and the accelerated medium is much different 
from the usual case (i.e., a fireball ejecta interacting with 
the static medium) . So the presence of neutron ejecta qual- 



itatively changes the early dynamical evolution of GRBs' 
remnant. 

In PD02 and B03, only the "fast" neutron component 
has been taken into account. In fact, the long, com- 
plex GRBs are more likely to be powered by the inter- 
action of shells with variant LFs, i.e., the internal shocks 
model (Paczynski & Xu 1994; Rees & Meszaros 1994). 
The best fit to the multi-wavelength afterglows implies 
that a significant fraction of the initial kinetic energy has 
been converted into internal energy (Panaitescu & Ku- 
mar 2002), which requires that the difference in veloci- 
ties between the shells is significant (i.e., the correspond- 
ing LFs satisfy rjf ^> rj s ) and their masses are compara- 
ble (m f w m s ) (Piran 1999, hereafter P99) [Throughout 
this work, the subscript f, s represent the fast/slow shells 
respectively; n, p represent the neutron/proton compo- 
nent respectively]. Typically, rj s is in order of tens, and 
rj{ is in order of hundreds. Thus, for the slow shell, 
the n, p component coast with r n s = r p s ~ rj s since 

r? s < Jfr ~ 3.9xl0 2 L^ 4 ro;7 /4 [(l+O/2]~ 1/4 (BM00), where 
L being the total luminosity of the ejecta, ro being the ra- 
dius of the central engine, £ being the ratio of the neutrons 
to the protons contained in the shells. In this Letter, we 
adopt the convention Q x = Q/i0 x for expressing the phys- 
ical parameters, using cgs units. For the fast shell, the n, 
p components move with different LFs: r nj f < r p f since 
generally r){ > rj^ (BM00). In the internal shocks phase, 
the fast ion shell catches the slower but earlier one at a 
radius R- m t ~ 10 13 cm and forms a new one moving with 
a LF r m ~ several hundred (hereafter the formed "new" 
shell is named as the "i-shell"). At i?i„t, the (3— decay of 
the neutron component is far from important unless the 
typical variability timescale is longer than 0.1 s (See Rossi, 
Beloborodov & Rees 2004 for more information). 

As the duration of GRB is long enough, the much later 
ejected i-shell catches up the earlier slow neutron shell 
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(hereafter, the "n-shell") at i? cat « 2Tl s c5T/ \l+z), where 
ST being the ejection time-lag of the earlier slow n-shcll 
and the later i-shell; On the other hand, the (3— decay ra- 
dius of slow neutrons reads Rp iS 2.6 x 10 13 r niS cm, at 
which slow neutrons have decayed significantly. As long as 
i? ca t is comparable with Rp, s , i.e., ST > 14(1 + z)T~], 5 , 
the decayed products of earlier but slower ejected neutron 
shells will be swept orderly by the later ejected ion shells 
in a range ~ i?; nt — Rp, s - The possible emission signature 
powered by that interaction is of our interest in this Letter. 
Otherwise, for GRBs much shorter than 14(1 + z)T~l 1 5 s, 
the last i-shell crosses the first slow n-shell at a radius 
<C Rp, s , where the (3— decay is unimportant and there is 
no observable signature. 

2. THE SIMPLEST MODEL 

Even for a neutron-free internal shocks model, the prac- 
tical LF distribution of shells involved is far from clear so 
far, let alone say the neutron-rich one. Here we use an ap- 
proximate model proposed by Guetta, Spada & Waxman 
(2001), in which the LFs of i-shells are drawn from a bi- 
modal distribution, r c j = r]f or T c j = ?7 S , with equal proba- 
bility. That simple model is favored by its relative high en- 
ergy conversion efficiency and peak energy. In this Letter, 
mainly for simplicity, we make the following assumptions: 
(i) £ = 1, and the ratios between the mass of the fast neu- 
trons, slow neutrons and i-shells are 1:1:2. (ii) Each fast 
(slow) neutron/ion shell moves with a same LF T ni f ~ 300 
(r n>s - 30) and r P)f - 1000 (r PiS ~ 30) respectively. After 
the merger of a pair of fast /slow ion shells, the resulting 
i-shell moves with r m ~ 200. (iii) The energetic ejecta 
are expanding into the low density interstellar medium 
(ISM), which is the case for most GRBs (Panaitescu & 
Kumar 2002). (iv) The duration of the GRB is longer 
than 30(1 + z)s since the constraint ST ~ 14(1 + z^"^ s 
should be satisfied. 

The (3— decay radius of fast neutrons Rp.f w 8 x 
10 15 r ni f 2 5 cm is much larger than J?/3 lS . We ignore them 
since for R < -R/3 >s , there is only a small fraction of fast 
neutrons decayed. The ISM has been swept by the early 
and fast ions as well as the decayed products of fast neu- 
trons, so the decay products of slow neutrons move freely. 

Assuming that at radius R, the (k — l)th i-shell crosses 
the jth slow n-shell (j, k 3> 1). At radius ~ R + AR, 
the fcth i-shell will catch up the jth slow n-shell, where 
AR w 2r 2 S A, A w cSt/(l + z), St ~ 0.01(l + z)s being the 
typical variability timescale of GRBs. For the jth n-shcll, 
the products of slow neutrons decayed at radius < R have 
been carried away by i-shells ejected early than the fcth 
one. Hence the mass swept by kth i-shell is AM ~ M r^ R i 

where M n = M° exp(— 7^-), being the the initial rest 

mass of the n-shell. So the averaged comoving number 
density of proton (due to (3— decay) swept by fcth i-shell 
can be estimated by 
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Notice that in front of the much later ejected i-shells, 
there are hundreds of decaying n-shells. With assumptions 
made before, the whole process of each i-shell interacting 
with these decaying n-shells is rather similar. For conve- 
nience, in our following treatment, the discrete interaction 



of each i-shell with these decaying n-shells has been sim- 
plified as an i-shell sweeping a moving proton trail (the 
LF of which is r n s ) with a number density n continually. 
The dynamics of that interaction can be described by the 
shock model as follows: Now the trail is moving with a 
LF r njS ~ tens, the generated thermal energy (comoving 
frame) in the shock front satisfies dU = (7 re i — l)dm c 2 
(e.g., B03) rather than dU = (7 — l)dm c 2 ; where 7 ro i = 
7r n , s (l — P-yPv^s) being the LF of the decelerating i-shell 
(the LF of which is 7) relative to the trail; (3 1 and /3r n „ 
being the corresponding velocity of 7 and r niS ; the mass 
swept by the i-shell reads 

dm = ATmm p T n . s R 2 ((3 7 - (3 Tne )d R ~ ^-dR. (2) 

After some simple algebra, the energy conservation of the 
system (the decelerating i-shell and the swept slow neutron 
trail) yields 
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where M; on being the rest mass of the i-shell; e being the 
radiation efficiency of the shock. For r n>s = 1, equation 

(3) reduces to the familiar form = _ 

(P99). 

In the downstream, the electrons have been heated by 
the shock. As usual (e.g., P99), we assume that the 
shocked electrons distribute as cx 7~ p for 7 C > 7 c , m , 
where p' ~ 2.3 being the typical power law index, 7 Cjm w 

e e^75f ^ E (7rci — 1); £ o and ee being the fractions of the 
shock energy given to electrons and the random mag- 
netic field at shock, respectively. The comoving down- 
stream magnetic field B w -\/327reB7rei(7rei — l)nm p c 2 . 
As shown in P99, there is a critical LF, above which 
synchrotron radiation is significant: 7^ = ( 1+ y™ T c ^ B 2 t , 

where t = i bs/(l + z) is determined by dR = 2j 2 cdt 
(t bs being the observer time, z <~ 1 being the red-shift of 
GRBs), <7t being the Thomson cross section, Y being the 
Compton parameter (e.g., Wei & Lu 1998, 2000; Sari & 

-l + A /l+4a:£ c /eB 



Esin 2001): Y 



where x being the radi- 



ation coefficient of electrons, so e = xe c . For 7 e , m > 7c, c 

p'-2 

(fast cooling), x = 1; For 7 Cim < 7^, x = (7c,m/7c,c)~ 
(slow cooling). 

Given a proper boundary condition, equation (3) can 
be easily solved. With resulting 7 , m , 7c, c and B, we 
can calculate the typical synchrotron radiation frequency 

v m = tttt^t— (e being the charge of electron), the cool- 

m 2(l+z)7rm c c v o e> n 



ing frequency v c 



7e,c7e-B 



and the self-absorption fre- 



2(l+z)-7rm c c 

quency f a (the detailed calculation on v a can be found 
in the appendix of Wu et al. (2003)), with which we 
can analytically calculate the synchrotron radiation flux 
at the fixed band. The numerical results for these typ- 
ical frequencies have been plotted in figure 1, where the 
parameters are taken as: M- lon = 2M° = 5.6 x 10 26 g, 
which is corresponding to a Luminosity L — 10 52 ergs s _1 
and St = 10~ 2 (1 + z)s, T m = 200; T n . s = 30, e c = 0.3, 
e B = 0.01, z = 1 (£> L = 2.2 x 10 28 cm). As shown in figure 
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1, at the early time, the electrons are usually in fast cooling 
case and u a is above the optical band. Lately, the electrons 
are in slow cooling case and v & drops to ~ a few x 10 13 Hz 
(Currently, the observed emission are coming from hun- 
dreds of i-shclls interacting with the slow neutrons trail. 
For the same reason, in our estimating the self-absorption 
frequency, the total number of electrons contributed has 
been assumed to be 300 times that of one trail. In fact, if 
just an i-shell interacting with the slow neutron trail has 
been taken into account, v a should be one order or more 
less than the value presented here). 

The synchrotron flux as a function of observer frequency 
can be approximated as follows: In the case of fast cooling 
(Just for z/ a > min{z/ c , z/ m }- For v a < mm{v c , z/ m }, please 
refer to Eqs. (4-5) of Zhang & Meszaros (2004)), 
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In the case of slow cooling 
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{V < Vc), 

(v c < V < v a ), 
(v a < V < v m ), 

{v > I/ m ). 

(4) 



{v < u m ), 
(v m < V < u a ), 
(u a < V < v c ), 
{v > v c ). 



(5) 

Approximately, F max = p 327r 2 e£ >2 jB, $ p , 

being a function of p' (for p' <~ 2.3, <I>p' ~ 0.6); N c be- 
ing the number of electrons involved in the emission, .Dl 
being the luminosity distance (Wijers & Galama 1999). 

The upcoming Satellite Swift carries three telescopes 1 : 
The Burst Alert Telescope (BAT), the X-ray Telescope 
(XRT), the Ultraviolet and Optical Telescope (UVOT). 
The energy range of XRT is 0.2-10 keV. UVOT covers 
170nm-650nm with six colors. So in Figure 2 we cal- 
culate the emission at the observer frequency z^ b s .i = 
1.0 x 10 15 Hz and v ohsa = 5kcV. 

The synchrotron radiation of an i-shell interacting with 
the trail as a function of the observer time has been shown 
in Figure 2. What we observed is the emission comes from 
a series (several thousand) of i-shells interacting with the 
trail of slow neutrons rather than just from one. In the 
simplest model proposed here, ~ 3s(l + z)/5t — 300 sam- 
ple light curves overlap, each with a St/(l + z) = 10~ 2 s 
delay. The resulting net observable fluxes at 10 15 Hz and 
5keV are about O.llJy and 1.4 x 10~ 9 ergs cm _2 s _1 respec- 
tively. The signature of 5keV emission is only marginal 
since the typical X-ray component during many GRBs is 
strong, too. So we do not discuss it any longer. The sen- 
sitivity of UVOT is about 19th magnitude in a 10-s expo- 
sure. The spacecraft's time-to-targct is about 20—70 sec- 
onds, in principle, part of, if not all the emission predicted 
here can be detected in the near future. However, we need 
to investigate whether some other emissions during prompt 
7— ray emission phase can dominate over the emissions 
predicted here or not. Recently, it is suggested that the 
relativistic e pair generated in the internal shocks phase 

1 http:/ /swift. gsfc.nasa.gov/docs/swift/ 



can power a UV flash ~ 13th magnitude, too (Fan & Wei 
2004). But the actual value of that emission is propor- 
tional to T 4 n [i5i/(1 + z)] 5 / 2 , and in that work the fiducial 
value for these parameters are taken as 300 and 0.1(1 + z)s 
respectively. If St ~ 0.01(l+z)s or smaller, as taken in this 
work, the UV flash predicted in Fan & Wei (2004) will be 
much weaker (<J 18th magnitude) unless T m is much larger 
than 300. Thus the UV emission predicted in this work 
can be detected independently. 

3. DISCUSSION 

In the standard baryonic internal shocks model for 
Gamma-ray bursts (Paczyhski & Xu 1994; Rees & 
Meszaros 1994), the prompt 7— ray emission are powered 
by the interaction between ion shells with variable LFs. 
If these shells are neutron rich, the neutron component 
coupled with the "fast ions" is accelerated to a larger 
bulk LF (~ several hundred), i.e., the fast neutron compo- 
nent, which is the focus of many publications (e.g., PD02, 
B03). However, the neutron component coupled with the 
"slow ions" can only be accelerated to a moderate bulk LF 
~ tens, for which the typical j3— decay radius ~ 10 15 cm is 
significantly smaller than that of the fast neutrons. As 
long as the duration of GRB is long enough, the earlier 
n-shells will be swept by the much later i-shells succes- 
sively till the radius ~ R/3, s , where the slow neutrons have 
decayed significantly. Consequently the interaction may 
power observable emission. 

In this Letter, with an extremely simplified model to 
describe the later but faster i-shclls interacting with the 
much early but slow, decaying n-shells, we have shown 
that there comes a <~12th magnitude UV flash at the later 
time of long GRBs. That emission is bright enough to be 
detected by the upcoming Satellite SWIFT in the near fu- 
ture. The emission predicted here is independent on the 
poorly known medium distribution around the progenitor 
of GRBs, as long as it is not very dense; it is also inde- 
pendent of some other physical processes, such as the e ± 
pair emission (Fan & Wei 2004), the possible reverse shock 
emission (As usual, if we assume the number density of 
interstellar medium (ISM) is about 1 cm~ 3 , Rp^ is much 
smaller than the deceleration radius ~ 10 17 cm, at which 
about half of the kinetic energy of the ejecta has been con- 
verted into the thermal energy of the shocked ISM, for a 
neutron-free ejecta) and so on. 

The resulting net light-curve in this work increases 
rapidly firstly. Then there comes a flat lasting ~ T90, the 
duration of corresponding GRB. Finally, the light-curve 
drops sharply, where the observed emission are mainly 
contributed by the "Equal-arriving surface" effect (e.g., 
Kumar & Panaitescu 2000; Fan, Wei & Wang 2004), if 
other emission such as the reverse shock emission has not 
been taken into account. It should be noted that the model 
taken here is a great simplification of the real situation, in 
which both the ion/neutron shells are moving with vari- 
able LFs, it is natural to expect that the corresponding 
emission is variable with time, too. However, to simu- 
late such complicate processes is far beyond the scope of 
this Letter. Nonetheless, the actual flux may not be far 
dimmer than our prediction since averagely speaking, the 
LFs of the slow n-shells and the i-shells are significantly 
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different, but their masses are comparable. So the inter- 
action of them can power a bright UV flash in the range 
- fO 13 - 10 15 cm. 
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Fig. 1. — The radiation frequencies involved in an i-shcll interacting with the trail of slow n-shells as a function of observer time. The 
physical parameters are taken as: Mj on = 2Mj = 5.6 X 10 26 g, which is corresponding to a Luminosity L = 10 52 ergs s" _1 and St = 10- 2 (l + z)s, 
r m = 200; r n , s = 30. e c = 0.3, <eb = 0.01, 2 = 1 (_D L = 2.2 X 10 28 cm). The solid line, dashed line and dotted line represent u c (the cooling 
frequency), v & (the synchrotron self-absorption frequency) and i/ m (the typical synchrotron frequency), respectively. The start point has been 
chosen to be R = 10 14 cm, the corresponding start time t b s ,beg = 0.04(1 + z)s. The end point has been chosen as 2.5 X 10 15 cm ~ 3Rp tS , the 
corresponding end time t b s .end — 3(1 + z)s. 
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Fig. 2. — The synchrotron radiation powered by an i-shcll interacting with the trail of slow n-shells as a function of observer time. The 
parameters are the same as those taken in Fig. 1. The solid line, dotted line represent u a ^ s = 10 15 Hz and 5 keV respectively. For clarity, 
the value of dotted line has been multiplied by a factor 100. Please note that what we observed is the emission powered by a series (several 
thousand) of i-shells interacting with the slow neutrons' trail rather than just by one. In the simplest model proposed here, about 300 
sample light curves overlap, each lags with St = 10~ 2 (1 + z)s. The resulting net observable fluxes at 10 15 Hz and 5keV are about O.llJy and 
1.4 X 10 _9 ergs cm~ 2 s _1 respectively. 



